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I , ABSTRACT 

Several massive black holes exhibit flux variability on time scales that correspond to 
source sizes of the order of few Schwarzschild radii. We survey the potential of near- 
infrared and X-ray polarimetry to constrain physical properties of such black hole 
systems, namely their spin and inclination. We have focused on a model where an 
orbiting hot spot is embedded in an accretion disk. A new method of searching for 
the time-lag between orthogonal polarization channels is developed and applied to an 
ensemble of hot spot models that samples a wide range of parameter space. We found 
that the hot spot model predicts signatures in polarized light which are in the range 
to be measured directly in the near future. However, our estimations are predicted 
^—Vj ' upon the assumption of a Keplerian velocity distribution inside the flow where the 

dominant part of the magnetic field is toroidal. 

We also found that if the right model of the accretion flow can be chosen for 
each source (e.g. on the basis of magnetohydrodynamics simulations) then the black 
hole spin and inclination can be constrained to a small two-dimensional area in the 
spin-inclination space. The results of the application of the method to the available 
near-infrared polarimetric data of Sgr A* is presented. It is shown that even with the 
currently available data the spin and inclination of Sgr A* can be constrained. Next 
generations of near-infrared and X-ray polarimeters should be able to exploit this tool. 
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1 INTRODUCTION 

The high-frequency emission from black-hole systems origi- 
nates mostly close to the black-hole horizon. In this central 
engine, huge amount of gravitational potential energy is re- 
leased from gas falling into the black hole. It is known that 
this emission can vary extremely rapidly (e.g. see Gaskell et 
al. 2007). This variable emission from compact objects has 
been interpreted by several authors as being related to the 
presence of orbiting spots in the inner parts of the accretion 
flows (Karas & Bao 1992; Karas 1999; Shahbaz 1999; Stella 
& Vietri 1999; Dovciak et al. 2004; Schnittman et al. 2006; 
Broderick & Loeb 2006, Zamaninasab et al. 2010). Here we 
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will focus on a new method based on the high-frequency 
time-resolved polarimetric observations which can constrain 
the intrinsic parameters of the sources. 

In the case of stellar-mass black hole candidates, the 
resemblance between the quasi-periodic oscillation (QPO) 
time scales and the dynamical frequencies of the black holes 
near their innermost stable circular orbit (ISCO) has in- 
spired the hot spot model (Stella & Vietri 1998, 1999). An 
orbiting blob inside the accretion flow will produce observ- 
able effects in both total and polarized flux. In the case of 
super-massive black holes, the situation is even more chal- 
lenging because the measured fluxes are much lower and 
time scales are much longer. Nonetheless, it is expected that 
similar physical mechanisms act in both categories of ac- 
creting black holes and, after appropriate rescaling, can be 
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Figure 1. Apparent images of an orbiting spot inside the ac- 
cretion flow around a Kerr black hole with the spin a = 0.5 
viewed by an observer inclined by 35°. The image on top shows 
the ray-tracing output for the total flux, while the middle and 
bottom show the flux in the 0° and 90° polarization channels, 
respectively. Here we have assumed that the projected spin axis 
is aligned toward the north on the observer's sky (8 = 0°) 



considered by similar techniques (Karas and Matt 2007). In 
Zamaninasab et al (2010), we have discussed in detail how 
flux magnification and changes in polarization angle (degree) 
follow specific patterns in this scenario, and we probed such 
a pattern in a sample of the near-infrared (NIR) flares of 
Sgr A*. 

Several mechanisms have been proposed for the cre- 
ation of the hot spots, including the random variations re- 
sulted from magnetic turbulence inside the magnetohydro- 
dynamics (MHD) flow (Balbus & Hawley 1991, Armitage 
& Reynolds 2003), vortices and flux tubes (Abramowicz et 
al. 1992), magnetic flares (Poutanen & Fabian 1999, Zycki 
2002) or star-disk interactions (Nayakshin et al. 2004, Dai 
et al. 2010). It is also shown that multi-component spot sce- 
narios are able to reproduce the overall behaviour of the 
observed power spectral densities (PSDs) and their (tran- 
sient) QPO features (Schnittman et al. 2006, Pechacek et 
al. 2008, Zamaninasab et al. 2010). 

In this paper, we will discuss a straightforward method 
which can be used for constraining the physical parame- 
ters of black holes (or other compact objects) which show 
variable polarized activity. The new method has several 
advantages in comparison with previous works (Meyer et 
al. 2006a, b, 2007, Hamaus et al. 2009, Zamaninasab et al. 
2008a,b, 2010). We must note that this method, in its core, 
assumes that variabilities in the flux are caused by orbit- 
ing spots. Normally, it has been reported in the literature, 
that such spots happen to be located at the ISCO. Here we 
will consider other possibilities for the spot's location and 
show how they will affect the predictions of the model (see 
Dovciak et al. 2004, Broderick & Loeb 2006 and Reid et al. 
2008 for detailed discussion about the spot model and ef- 
fect of the position of the spot on the observed signal). We 
must also stress that the results derived here are dependent 
on our assumptions about the Keplerian velocity of the flow 
and the magnetic field orientation. 

In section [5] we will describe how the spot model pre- 
dicts a time-lag between orthogonal polarization channels 
and study how these time-lags are characterized by the black 
hole mass, spin and inclination angle. A description about 
the available NIR data from the Galactic Centre black hole 
and how they can be used to constrain its spin and incli- 
nation is provided in section [3] We also briefly discuss a 
possible way to apply this method to other sources, like 
RE J1034+396. In section]!] we summarize the main results 
and draw our conclusions. 



2 TIME DELAY BETWEEN POLARIZATION 
CHANNELS: PREDICTIONS OF THE HOT 
SPOT MODEL 

MHD simulations have demonstrated that in the inner parts 
of the accreting plasmas, magnetic-field dissipation mod- 
erates the field intensity far below its equipartition value 
(Kowalenko & Melia 1999). Nevertheless, when the gas prop- 
agates and follows a Keplerian flow, a MHD dynamo can 
generate an increased (yet subequipartition) magnetic field 
that is dominated by its azimuthal component (Hawley, 
Gammie, & Balbus 1996). In this paper we examine fluc- 
tuations occurring in the inner parts of the black hole's ac- 
cretion flow, very close to the event horizon. Since our focus 
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Figure 2. Top: Total flux magnification and the corresponding 
polarization channels for a spot orbiting at the marginally stable 
orbit of a Schwarzschild black hole and viewed by an observer in- 
clined by 60° . Flux densities are all normalized to the maximum 
of the total flux value. Bottom: Cross-correlation between light 
curves of orthogonal channels: 0° & 90° (black) and 45° & 135° 
(red). The cross-correlation shows a clear delay between orthog- 
onal pairs of channels with the values of 5i and &2 ■ Here we have 
assumed that the projected spin axis is aligned toward the north 
on the observer's sky (8 = 0°) 



is mainly on radii smaller than 40 r g (where r g = is the 
gravitational radius) we have adopted a thin disk approxi- 
mation while magnetic field is dominated by its azimuthal 
component. Thin disk approximation was also chosen for 
reasons of practicability. It allows in a straight forward way 
(low computational effort) to repeat simulations for a wide 
range of free-parameters space. Adopting a relativistic treat- 
ment of the inner Keplerian region, while considering certain 
aspects of earlier hydrodynamical simulations, we end up 
with a model where the gas circulates at small radii twist- 
ing the magnetic field inside a relativistic Keplerian disk. To 
estimate the possible structure of the gas distribution within 
several Schwarzschild radii of the black hole, we followed the 
investigation done by Qian et al. (2009). 

For the radiation mechanism, we assumed that NIR and 
X-ray synchrotron photons originate from the relativistic 
electrons gyrating around magnetic-field lines inside the flow 
(see Liu et al. 2006, Eckart et al. 2004, 2006, 2008, Zaman- 
inasab et al. 2010). The energy distribution of electrons is 
approximated by a power-law function: 
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where N("/), 7 and y c are electron energy distribution func- 
tion, Lorentz factor of the electrons and its cut-off, respec- 



tively. Thus we can define both the ordinary and extraor- 
dinary emission coefficients for each point of the accretion 
disk (Pacholczyk 1970; see also Zamaninasab et al. 2010): 
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if 5/3 and K 2 /3 are the corresponding modified Bessel func- 
tions, v is the frequency measured in the co-moving frame, 
B is the magnetic-field strength, and 8 e is the angle between 
the direction of the magnetic field and the direction toward 
the co-moving observer: 
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In the above equation B is the magnetic field four-vector and 
p e represents the direction of the emitted photon momentum 
in the co-moving frame. 

Since the majority of the high-frequency photons orig- 
inate from plasma very close to the black hole, the curved 
structure of space-time must be taken into account. In our 
simulations, we have used Karas-Yaqoob (KY) code for the 
ray-tracing (see Dovciak et al. 2004 for details of the ray- 
tracing method). Light bending, aberration, Doppler boost- 
ing, and frame dragging affect both flux and polarization 
properties of the emission. The KY code is able to calculate 
all the relevant relativistic effects for the Kerr space-time. 
The observed flux density from the accretion disc at a cer- 
tain observed frequency, F Vo , can be computed as follows 
(Dovciak et. al. 2004b, 2008a,b): 
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(7) 



where I v is the intensity, D is the distance to the source, 
dA = rdrdtf) is an area element on the disk and G(r, (f>) is 
the transfer function which contains all the above mentioned 
relativistic effects. This function depends on the location of 
the emitter in the disc. 

The hot spot is modelled as an over-density of non- 
thermal electrons centred at a point orbiting at the Kep- 
lerian velocity with a Gaussian profile as measured in the 
co- moving frame: 

— (Ar) 2 

n = n exp[— ], (8) 

where Ar = |r — with ~r s p being the vector pointing 
to the centre of the spot. We have chosen the typical value 
size of one Schwarzschild radius for the size of the spots. For 
the underlying steady accretion flow, we assumed that both 
non-thermal electron density and the magnetic field strength 
scale as r~ . Since our analysis is focused on high-frequency 
regimes (NIR and X-ray), we have ignored radiative-transfer 
effects. 

Here we must note that some models propose strong 
deviations from Keplerian velocity for the accretion flow. 
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Figure 3. Time lag between the orthogonal polarization channels 5i(= <5o— 90), &2(= <545-i3s) and flux magnification /i, as functions of 
spin and inclination (from left to right, respectively). The top row corresponds to § = 1.0 while the bottom row shows the same plots 
for £ = 1.5. Here we have assumed that the projected spin axis is aligned toward the north on the observer's sky (9 = 0°). Values are 
presented in gravitational unit time scale (t g ). 



In the absence of an unambiguous theory for the accretion 
model, we have adopted the simple case of the Keplerian ve- 
locity distribution. While our choice of a thin accretion disk 
reduced the computational time needed for the ray-tracing 
simulations and made it possible to cover a wide range of free 
parameters in a reasonable time, it can produce large dis- 
crepancies at high inclinations. For this reason, we have lim- 
ited ourselves to inclinations less than 80°. The dimension- 
ality of the spot is also important for polarized ray-tracing 
along strongly lensed rays. In this case, many different mag- 
netic field orientations could be encountered (as seen by an 
observer who is parallel propagating his orientation along 
the rays). For high spin values this can produce substantial 
depolarization. Our comparison with three-dimensional sim- 
ulations with the same setup, showed that for inclinations 
less than 75° and black hole spin below 0.9 the deviations 
of the models are less than 15% while above this values it 
can reach as high as 40%. 

Figure [T] shows the images of the assumed accretion 
flow at an ideal detector as well as polarimetric images at 
0° and 90° polarization channels for a spinning black hole 



(a — 0.5) viewed by an observer inclined by 35°. a and /3 are 
the projections of the impact parameter of the emitted pho- 
tons onto the sky. Here we have assumed that the projected 
spin axis of the black hole coincides with the north on the 
observer's sky. We have used these images to extract light 
curves of the total flux and the corresponding polarization 
channels. Figure[2]shows the behaviour of flux light curves in 
4 orthogonal polarization channels (namely 0° , 45° , 90° and 
135° channels) for an orbiting blob located at the ISCO of 
a Schwarzschild black hole, viewed by an observer inclined 
by 60°. One can see clearly that each channel reaches its 
maximum at a different time. The reason for this behaviour 
is the strong Doppler beaming as well as the light focusing. 
Furthermore, the rotation of the polarization plane along 
the photon trajectory plays a role. The latter effect is par- 
ticularly strong for small radii of the spot orbit (Dovciak 
et al. 2008a,b). As Fig. [2] shows, the time-lag between each 
orthogonal pair of polarization channels can be measured by 
cross-correlating the corresponding light curves. One can de- 
fine another parameter called magnification factor (n) which 
provides useful information about the properties of a light 
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Figure 4. Normalized probability of creation of spots inside the 
accretion disk as a function of distance from the black hole (£). 
Different curves correspond to different models: (A) spots always 
created very close to the ISCO, (B) dominant creation happens at 
ISCO and gradually decreases with increasing radius, (C) maxi- 
mum probability happens at 2 X ISCO and (D) equal probability 
everywhere inside the disk. 
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Figure 5. Probability densities for a given inclination and black 
hole spin marginalised over 6 for the hypothetical values of <5i = 
(4±0.5) t g , &2 = (3±0.5) t g and = (4±1). This example shows 
how high accuracy measurements of time-lags and magnification 
can result in tight constraints on a and i. 



curve. Here, the magnification factor refers to the ratio of 
the maximum flux density to its quiescent value. This factor 
is a measure of how much the flux is magnified by lensing 
and boosting effects. 

Up to here, we assumed that the projected spin axis of 
the black hole is aligned toward the north on the observer's 
sky. In general, the measured time-delay between two or- 
thogonal polarization channels depends on the orientation 
of the projected spin axis with respect to the north (9, van- 
ishing at the north and increasing toward the east) in our 
chosen reference coordinates. In order to reduce this depen- 
dency, one can measure simultaneously time delays between 
two pairs of orthogonal channels which are rotated by 45° 
with relative to one another. Since changing the position an- 
gle corresponds simply to a coordinate rotation, the images 
were originally computed only for a single position angle, 
namely 9 — 0°. The time-lag associated with each desired 
value of 6 were extracted from these images by applying the 
appropriate Mueller matrices. 

The main question which rises in this step is how these 
time lags (5i = <5o-90 and 82 = 545-135) and magnification 
factor (fi) are sensitive to the changes of the global parame- 
ters of our model. The main free parameters are the spin of 
the black hole (0 ^ a 1 in geometric dimensional units), 
inclination of the observer (0° ^ i ^ 90°, when the edge-on 
view corresponds to i = 90°), orientation of the projected 
spin axis on the sky (—90° ^ 9 ^ 90°, vanishing at the 
north and increasing toward the east), and the distance of 
the blob from the black hole (controlled by the parameter 



f x ISCO). Figure [3] shows Si, 82 and fi as func- 



tions of the spin and inclination for a black hole of the mass 
M in gravitational units. The functions have been presented 
for two distances of the spot from the black hole (£ = 1.0 
and £ = 1.5) while 9 = 0°. We have calculated these func- 
tions for a wide range of acceptable values of this parameter 
(£ = 0.5 — 3.0; with a step of 0.1). One can see that these 
three different functions behave differently in the a — i space. 
This allows us to use these functions as tools for constraining 
spin, inclination and orientation of the accreting systems. 

It is clear from Fig. [3] that our results are sensitive to 
the exact location of the orbiting blob with respect to the 
black hole (see predictions for £ = 1.0 and £ = 1.5 in Fig.[3j). 
Unfortunately, there is no direct way to constrain that pa- 
rameter. A lot of effort have been devoted to understand 
where the radiation and the stress edges of an accretion 
disk are located and how they depend on the different as- 
sumptions of the models (Beckwith et al. 2008, Shafee et al. 
2008, Hilburn et al. 2010, Noble et al. 2009, 2010). Since it 
is not yet clear how and where spots can be created inside 
an accretion disk, we have assumed four different types of 
models for the spatial distribution of the spots as a function 
of distance from the black hole (Fig. [4}. We have considered 
models in which the main creation of spots happens only 
around the ISCO (model A), distributed over the disk with 
the maximum at the ISCO and gradually decreases with 
increasing the radius (model B), distributed over the disk 
with the maximum at the 2x ISCO (model C, see Shafee 
et al. 2008) and an equal distribution for the probability of 
a spot being located anywhere inside the disk (model D). 
We have used the product of the magnetic stress and energy 
dissipation functions from MHD simulations (e.g. Hilburn 
et al. 2009 and Shafee et al. 2008) in order to define such 



6 Zamaninasab et al. 





0.0 0.2 0.4 0.6 0.8 



20 40 60 80 



i [deg] 




Figure 6. Probability densities for a given spin p(a), inclination p(i) and orientation of the system on the sky p(0), marginalised over 
other parameters (left to right, respectively). The assumed values of Si, 62 and fi are the same as in Fig. [5] The dark and light shaded 
areas denote the lcr and 2cr confidence limits, as defined in the text, respectively. 



radial distributions. We have adopted a Bayesian parameter 
estimation method used by Broderick et al. 2009. Here, we 
repeat the basics of the method in order to keep the pa- 
per self-contained (see Broderick et al. 2009 for more details 
about the method). 

The constraints upon a, i and 8 we report here, suffer 
from the uncertainty introduced by the assumption that our 
model is appropriate for the observations. We compute the 
probability following a Bayesian scheme, which states that 
a given set of model parameters (a,i,8) are correct for the 
measured delays and magnifications. 

Given a particular set of model parameters (a,i,#, £), 
the probability of measuring an observable Xj (where Xj G 
[8ij, S2j, Pi]), assuming Gaussian observational errors, is 



p(a,i,6,£\{xj}) 



1 

■ exp 



2-irAxi 



[xj - x{a,i,6,i)f 
2Ax* 



dxi 



(9) 



Therefore, the probability of simultaneously observing a 
measured set of [5ij , , fJ>j] for a certain model is 



P{{5i j ,&2j,^j}\a,i,9,Cl = 
P{5ij\a,i,e^) x P(S 2j \a,i, 



>,0 xP(^\a,i,6,0- 



(10) 



Contrary to eq. {TO}, which gives the probability that the ob- 
served quantities come from a certain model, choosing an ap- 
propriate set of priors on a, i, 8 and £ allows us to compute, 
via Bayes' theorem, the probability density of a set of model 
parameters given the observed values, i.e. p(a,i,8,£\{xi}). 
As such, we now turn to the problem of choosing these pri- 
ors. 

Assuming the isotropic probability of the black hole's 
spin orientation, a choice for the prior upon i and 8 can 
be made so that p(i,9) — sin(i) (Broderick et al. 2009). 
The prior on a is chosen to be uniform, i.e. p(a) — 1 since 
we lack a complete theoretical knowledge about the spin 
evolution of super-massive black holes. Finally, the prior on £ 
is set according to the models in Figure [4] Therefore, Bayes' 
theorem gives 



(11) 



P({x ] }\a,i,8,Qp(a)p(i,8)p(0 
da di dd d£ P({xj}\a, i, 9, £)p(a)p(i, 8)p(£) 

P({x J }\a,i,8,£)sin(i) 

da di dd d£ P({xj}\a, i, 8, £) sin(i) 



As it can be seen, this is a probability density in a 
four-dimensional parameter space, which makes it difficult 
to be visualized directly. Therefore, we construct a variety 
of marginalised probabilities for presentation and analysis. 
This is done by forming a pair of two-dimensional p(a, i) 
function with the probability densities marginalised over 8 
and £: 



p{a,i) 



dd dtp(a,i,8^\{x 3 }) 



(12) 



This is plotted in Fig. [5] for a hypothetical example of a 
source with the measured delays of Si = (4±0.5) t g (where t g 
is the gravitational time) and 82 = (3 ±0.5) t g while the flux 
is magnified by a factor of /i = (4± 1). In an observationally- 
constrained astrophysical case, it is difficult to determine the 
values of the time delays and magnification so precisely. As 
one can see, the model allows to exclude a wide range of 
possible a-i combinations and well constrains the possible 
range of inclination and black hole's spin. In order to dis- 
tinguish the probability distribution of each parameter, the 
marginalised one-dimensional probability densities are also 
constructed: 



p{a) 



dd di d£p(a, i, 8, £,\{xj}) 



(13) 
(14) 
(15) 



These are shown in Fig.[5]for the same hypothetical values of 
Si, 82 and /1 as in Fig.(6] For both Figs.[5]and[6]we have used 



p(i) — J dd da d£p(a,i,8,£\{xj}) 
p{8) — I di da d£p(a,i,8,t;\{xj}) . 



the commonly used assumption that the spot is located at 
the ISCO (model A for p(£)). It is evident that the allowed 
parameter space is highly non-Gaussian. As a result, the 
extraction of values and their associated uncertainties for the 
fitting parameters must be done carefully. These values will 
be highly correlated in all occurrences, and the systematic 
uncertainties due to assuming a particular accretion model 
will dominate the errors (and thus will not be reviewed again 
here) . 

For determining the la and 2a error intervals, we used 
an additive probability as defined in Broderick et al. 2009: 

P(>p) = / p(x)dx (16) 

where x is the parameter of the total probability distribu- 
tion. The above definition of P(> p) gives the probability 
associated with the parameter space region that has a prob- 
ability density above p. We set the 1-a and 2a contours to 
indicate the regions of p for which P(> p) = 0.683 and 
P(> p) = 0.954, respectively. These follow the normal defi- 
nition of la and 2a errors, while in our case the errors are 
strongly non-Gaussian. In Fig. the la and 2a confidence 
intervals are indicated by the dark and light shaded regions. 

In the next section, we will discuss the implications of 
our model for the observed variable NIR polarized emission 
from Sagittarius A*. 



3 OBSERVATIONS 

3.1 Near-Infrared Flares of Sagittarius A* 

At the centre of the Milky Way stellar motions allow us to 
firmly associate Sagittarius A* (Sgr A*) with a ~ 4x 10 6 M 
super-massive black hole (Eckart & Genzel 1996, Genzel et 
al. 1997, 2000, Ghez et al. 1998, 2000, 2003, 2005a,2008, 
Eckart et al. 2002, Schodel et al. 2002, 2003, 2009, Eisen- 
hauer et al. 2003, 2005, Gillessen et al. 2009). 

Recent radio, near-infrared and X-ray observations have 
detected variable and polarized emission and give detailed 
insight into the physical emission mechanisms at work in 
Sgr A* (e.g. Baganoff et al. 2001, 2002, 2003, Eckart et 
al. 2002, 2004, 2006, 2008, 2009, Porquet et al. 2003, 2009, 
Goldwurm et al. 2003, Genzel et al. 2003, Reid et al. 2004, 
Ghez et al. 2004ab, Eisenhauer et al. 2005, Belanger et al. 
2006, Hornstein et al. 2007, Yusef-Zadeh et al. 2006ab, 2007, 
2008, 2009, Marrone et al. 2008, Sabha et al. 2010). Sgr A* 
- in terms of Eddington luminosity - is the faintest super- 
massive black hole known. However, due to its proximity 
it is bright enough to be studied in great detail. With the 
possible exception of the closest galaxies, no extragalactic 
super-massive black hole with a similar feeble Eddington 
rate would be observable. 

The NIR data we investigate here were taken in the 
K s band with the NIR camera CONICA and the adaptive 
optics (AO) module NAOS (NACO) at the ESO VLT unit 
telescope 4 (YEPUN) on Paranal, Chik-Q on 13 June 2004 



1 Based on observations at the Very Large Telescope (VLT) of 
the European Southern Observatory (ESO) on Paranal in Chile; 
Programs:075.B-0093 and 271.B-5019(A). 
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Figure 7. The NIR flare event of Sgr A* observed on 13 June 
2004. Four different polarized channels are indicated with dif- 
ferent colours (top). In order to derive the time-lags a cross- 
correlation analysis has been performed. Maximum of two cross- 
correlation curves have been indicated with two different meth- 
ods (middle and bottom). The average values of <5i = (4.20 ± 
1.00) min, 6 2 = (2.12 ± 1.00) min and p, = (2.98 ± 1.00) derived 
from both analysis. 



(start of observation: 07:20:02 UT). For determining the lin- 
ear polarization characteristics of a time-varying source, a 
Wollaston prism, installed on the NACO NIR camera, per- 
mits simultaneous measurements of two orthogonal direc- 
tions of the electric field vector. It is combined with a rotary 
half-wave plate to allow for rapid alternating measurements 
of the electric field vector at different angles. Details of the 
data reduction and flux calibration is described in Eckart et 
al. (2006a), and Zamaninasab et al. (2010). 

Figure [7] shows light curves of a polarized flare from 
Sgr A* observed on 13 June 2004 in four polarization chan- 
nels. There is a clear time-lag visible between orthogonal 
channels which is bigger than the sampling rate of the ob- 
servations. Furthermore, the light curve is continuous and 
there is no gap according to the usual sky background mea- 
surements at the time of the flare. Each channel shows an 
approximate Gaussian behaviour while clear time-lags be- 
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Figure 8. Same as Fig.[5]for the time-lag and magnification values derived from the observed flare of Sgr A* on 13 June 2004 (Fig.[7Jl. 
The results are shown for different models of p(§): model A (top left), model B (top right), model C (bottom left) and model D (bottom 
right). 



tween maxima (specially between orthogonal pairs) are vis- 
ible. 

We have calculated the time lag between orthogo- 
nal pair of channels for the original (Fig. [7J middle) and 
first-order interpolated light curves (Fig. [7] bottom). For 
the cross-correlation function we have followed Alexander 
(1997). Results of both methods are well in agreement 
with each other. This leads to the estimated values of 
the time delays between orthogonal polarization channels 
5i = (4.20 ± 1.00) min.fo = (2.12 ± 1.00) min and a mag- 



nification value of p, = (2.98 ± 1.00). The errors in these 
estimations derived from the standard deviation of \ 2 fit- 
ting of Gaussian functions to the observed flux and cross- 
correlation functions. We must note that the sampling rate 
of this observation was ~ 2 min, so the level of confidence 
specially for the 82 value is not high. In order to improve 
these estimates we need high SNR images observed in po- 
larimetry with higher time resolutions. 

Although the time resolution of the available data is 
not better than 1.5 minutes, it is good enough to rule out 
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Figure 9. Same as Fig.[6]for the time-lag and magnification values derived from the observed flare of Sgr A* on 13 June 2004 (Fig. [J}. 
p(a), p(i) and p(8) are marginalised over all other parameters assuming model A for p(£). 
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Figure 10. Same as Fig. [9] assuming model B for p(£). 
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some parts of the a — i plane. Fig. [8] shows two-dimensional 
marginalised probability densities p(a, i) for the measured 
values of magnification and time delays. The Four different 
images correspond to different assumptions about the spatial 
distribution of the spots (p(£)). We see that while the most 
popular assumption that spots are located around the fSCO 



50° 



the more 



(model A) results in a = 0.5_ ' 4 and 
reasonable assumption which allows the spots to be created 
over a wider range of locations increases the estimated values 
for both inclination and black hole's spin. 

Figs. 151 1121 show the fully marginalised probability dis- 
tributions for different assumptions about p(£), respectively. 
However, we must remind the reader that these parameter 
estimates are strongly correlated, and thus they should be 
used with caution. Here we again use the additive probabil- 
ity, P(> p) to define the la and 2a intervals, shown in the 
panels of Figs. 1 91121 as the dark and light shaded regions, 
respectively. 

From the left panel of Fig. [9] the most likely spin value 
is a = 

°- 45 ai5 2 ( 1<T confidence). While we can rule out very 
low spins (a ^ 0.1) at the 2a level for our particular choice 
of the spin is otherwise weakly constrained. Left panels 
of Figs. 1911121 show that this estimation is highly dependant 



on the choice of p(£) and has a tendency toward higher spin 
values for models B-D. 

In contrast, the inclination is robustly limited towards 
moderate and high angles. The most likely inclination value 
is i = 55°^2o° (1°" confidence). It can be seen clearly in the 
central panels of Figs.lWT2"lthat face-on geometries (i ^ 30°) 
are convincingly ruled out. However, the general trend of all 
the considered models is toward highly inclined systems. 

Interestingly, constraining the position angle shows that 
all the assumptions result in a more or less similar val- 
ues (6 — 50^qo ) as can be seen in the right panels of 

Figs, ma 

Our limits upon the black hole's spin are in good agree- 
ment with previous estimates by Huang et al. 2009, Dexter 
et al. 2010 and Moscibrodzka et al. 2009 while it contra- 
dicts the low spin values favoured by Broderick et al. 2009. 
However, we remind the reader again that the constraints 
on spin are not strong enough to rule out any possibility. 

The limits upon the inclination are also in quite good 
agreement with previous efforts (Meyer et al. 2006, Markoff 
et al. 2007, Broderick et al. 2009, Falanga et al. 2009, Huang 
et al. 2009, Moscibrodzka et al. 2009 and Dexter et al. 2010). 
They all favour high inclinations using different approaches 
to the problem. Falanga et al. 2009, for example, find an in- 
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Figure 11. Same as Fig. [9] assuming model C for p(£). 
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Figure 12. Same as Fig. l9l assuming model D for p(£). 



clination of about 77° ± 10° (similar to the results by Meyer 
et al. 2006) by fitting Sgr A*'s flares with a model of hy- 
drodynamic instabilities. On the other hand, Markoff et al. 
2007 fit the long-wavelength observations with a different 
approach, a hydrodynamic jet. They also favour large incli- 
nations (i 75°). 

Unlike the inclination, our position angle estimate dis- 
agrees significantly with many previous efforts except Meyer 
et al. 2006 and Markoff et al. 2007. Our estimated position 
angle is marginally consistent with the second solution of 
Broderick et al. 2009 at the la level. 

However, we have noticed that the our estimated 6 co- 
incides very well with the value derived by Muzic et al. 
2010 where they used the structure of two bow shock stars 
(namely X3 and X7) in the vicinity of Sgr A* to trace a 
possible outflow scenario (Fig. I13|) . Our best bet (6 = 50°) 
exactly passes the positions of X3, X7 and the mini-cavity as 
shown in Fig. 1131 It is proposed that the whole mini-cavity 
structure can be caused by wind from the direction of Sgr A* 
implied by traces of hot gas and plasma shocks observed in 
this region (Yusef-Zadeh et al. 1990, Yusef-Zadeh & Melia 
1992, Yusef-Zadeh & Wardle 1993 and Muzic et al. 2007). 
At the same time, the results leave a lot of space for future 
improvements. 



3.2 The case of RE J1034+396 



The main advantage of this method is its applicability for a 
wide range of compact sources which show polarized variable 
emission in NIR and X-ray wavelengths. Since the current 
technical difficulties do not allow for polarimetric observa- 
tions with time resolutions less than ~ 30 seconds (at least as 
far as it concerns VLT), the method is currently applicable 
to sources with masses greater than ~ 1O 6 M0. This makes 
objects like RE J1034+396 perfect candidates for testing 
predictions of the hot spot model. The recent unambiguous 
discovery of one hour periodicity in the X-ray emission light 
curve of this source is interpreted to be related to the ISCO 
frequency of this super-massive black hole (Gierlifiski et al. 
2008) . Unfortunately, there is no polarimetric observation of 
this source available so far. The mass of the object is also 
not well determined (the estimates vary from 6.3 x 10 5 to 
3.6 x 1O 7 M0). This source looks to be a promising target 
for future X-ray polarimetry if its mass estimation can be 
improved (using reverberation mapping for example). 
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Figure 13. The estimated orientations of the Sgr A*'s spin axis 
in our sky from all the models of p(£) (6 = 50^jq) are m ver y 
well agreement with the study of Muzic et al. 2010 which relates 
the alignment of two bow shock sources (X3 & X7, marked with 
open and filled stars respectively) and mini-cavity as a possible 
hint for a wind from Sgr A*. 

4 SUMMARY AND CONCLUSION 

We have investigated the general relativistic effects on the 
polarization properties of the emission from an orbiting spot 
inside the accretion flow around a black hole. We have de- 
scribed a method in which observing the time delay between 
the polarization channels in high-frequency regimes can lead 
to a constraint on the spin and inclination of the black hole. 
The method may also be applied to the X-ray domain for 
sources which are bright enough in polarized X-ray flux. We 
applied the method to several hot spot models that sample 
a wide range of parameter space. We have discussed how 
the probability for the spatial distribution of the spots can 
affect our results. This shows the crucial need for reliable 
2 and 3-dimensional GRMHD simulations of accretion disks 
concerning the exact location of the plunging region and the 
stress edge. The results of the application of the method to 
the available NIR polarimetric data of Sgr A* is also pre- 
sented. The results for this source are consistent with previ- 
ous findings of an inclined accretion disk around a spinning 
black hole (a ^0.4). While with current VLT facilities the 
rate at which high S /N images in polarized flux density are 
taken can be increased significantly by a new method of ob- 
servation, the next generations of NIR and X-ray polarime- 
ters should be able to exploit this method as a powerful 
probe of the metric around black holes. 
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